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1. Abstract 



We study the stochastic quantization of the system with first class constraints 
in phase space. Though the Langevin equations of the canonical variables are 
defined without ordinary gauge fixing procedure, gauge fixing conditions are au- 
tomatically selected and introduced by imposing stochastic consistency conditions 
upon the first class constraints. Then the equilibrium solution of the Fokker-Planck 
equation is identical with corresponding path integral distribution. 
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2. §1. INTRODUCTION 

The stochastic quantization method, which has been proposed by Parisi and 
Wu [1) in 1981, is known as the third quantization method 121 following the canon- 
ical and path integral formalism. In the scheme fields (f>(x) have an additional 
coordinate t; 



coordinate which is called the fictitious time. The development of the fields with 
the lapse of the fictitious time is described by a Langevin equation 



where L is the Lagrangian of a system and dW is the Wiener process whose corre- 
lations are defined as 



The equations enable us to calculate correlation functions of the fields, which tend 
to corresponding Green functions in the equilibrium limit t — > oo. 

The correlation functions may be given in another form 



where P is a probability distribution. A Fokker-Planck equation, which P satis- 
fies, is obtained by calculating, with the Langevin equation, the correlation func- 
tion of the fictitious time derivative of any functional f((f>) and comparing it with 



</>(x) -> <j>(x,t), 




(dW(x,t)) = 0, 



(dW(x, t)dW(x', t)) = 25{x - x')dt. 
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J D(f)f ■ dP/dt. The Fokker-Planck equation thus becomes 

dP__d_rdL _cM 

The equation has in the equilibrium limit dPj dt = a solution 

P = exp(-L), 

which is the same as corresponding path integral distribution. 

It is one of the large merits of the stochastic quantization method that it enables 

[3] 

us to quantize gauge fields without gauge fixing conditions. Hence if we extend the 
stochastic quantization method, which is ordinary defined in configuration space, 
to phase space, we may shed light on the quantization of the systems with first 
class constraints (FCC). Until now Namiki et al. M and we M have studied how to 
handle constraints by the stochastic quantization method in configuration space 
(the SQM-cs), showing that it is much useful to impose the consistency condition 
with the lapse of the fictitious time, the stochastic consistency condition, upon 
the constraints. On the other hand, some proposals for the stochastic quantization 
method in phase space (the SQM-ps) have been madef 61 0hba l71 has tried to quantize 
by the SQM-ps, in particular, the system with FCC where gauge fixing conditions 
are added. It is accordingly the system with second class constraints that he 
quantized by the SQM-ps. 

In this paper we quantize the system with FCC by the SQM-ps without usual 
gauge fixing procedure and show that gauge fixing conditions are automatically 
introduced through the stochastic consistency conditions upon the FCC. Since the 
prescription to impose them is straightforward and has little dependence on the 
form of the FCC, it will be useful when we quantize complicatedly constrained 
fields and may bring us the other knowledge about gauge theory than the other 
quantization methods have offered. 

This paper is organized as follows. In §2, taking notice of the relation between 
the SQM-cs and the SQM-ps and, in particular, the Lee- Yang term, we define the 
strict form of the Langevin equations with a constant kernel, of canonical variables. 
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In §3, we show how to apply the SQM-ps to the system with FCC. We impose the 
stochastic consistency conditions upon the FCC and the determinant factor which 
appears with the gauge fixing procedure consequently becomes included in the 
solution of the Fokker-Planck equation in the equilibrium limit. In §4, we define 
new fields and divide them into constraint variables and independent ones to clarify 
what gauge the stochastic consistency conditions accord to. §5 is devoted to the 
conclusion. Besides them, we give in the appendix how to apply our formulation 
to the system with second class constraints (FCC and gauge fixing conditions). 

3. §2. THE LANGEVIN EQUATION IN PHASE SPACE 

In this section we consider the system with the Lagrangian in Euclidean space- 
time 

L = \G AB {q)d ll q A d ll q B , (2.1) 
A,B = l,-.-,N, 
(j, — 0, • • • , n. 

For such a system as includes a field- dependent metric, by using Ito's calculation 
rule, the Langevin equation in configuration space is 



dq A (t) = -G AB ^dt + ^-(y/GG^dt + E A a dW\ (2.2) 



PL _i a_ 

a? + 

where G AB is the inverse of the metric Gab and E A a is the vielbein defined with 
its inverse E A a as 

E A a E B b 5 ab = G AB , 
E/EsSab = Gab, 

rpA rp a jpA jp b r& 

E a E B = d B> E a h A = d a , 

a,b — 1, • • • , N. : reference frame indices 

The additional term of the Langevin equation (2.2) is necessary for general co- 
ordinate transformation (GCT) covariance l91 and, in perturbation, cancellation of 
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divergent terms.' 51 Moreover we can observe it induces the Lee- Yang term' 101 in equi- 
librium Fokker-Planck distribution!" 1 The Fokker-Planck equation in this case is 

dP d AB f dL 1 dVG d i 
Ot dq A \dq B dq B dq B J ' 

G = det(G AB ), 
the equilibrium solution of which is 

P = exp(-L + i/nG). 
The correlation function is thereby written as 

(•••) = j Dq---exp(-L + hnG), 

which is consistent with the path integral formalism. 

Let us shift to phase space. First we define the conjugate momentum p A of q 
and the Hamiltonian H as 



A 



PA = +i d ( doq A} = iG ABd q B , (2.3) 



H = +i PA d q A + L = \G AB p A p B + \G A Bdiq A diq B . (2.4) 
The equations of motion in phase space are the Hamilton equations 



dH 

-id p A - ^-j = 0, (2.5a) 



a dH 



+id q A - — = 0, (2.56) 
dpA 

in place of the Euler-Lagrange equation in configuration space. It is natural to 
adopt the LHS of the Hamilton equations (2.5) as naive drift terms of the Langevin 
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equations in the phase space. Moreover since the metric of the manifold spanned 
by the canonical variables q and p is 



( 



G AB {q) 
V G AB (q) 



determinant of which is unity, the GCT covariant Langevin equations are defined 
in our formulation as 

a tt a/~<AB 
dq A = G AB (-td 0PB ~ TTB)dt + -TTB-dt + E A a dW a , 

dq dq (2.6) 

r)J-f 

dp A = CG AB (+td q B - —)dt + VCE^dV a , 

dp B 

where C is the constant kernel determined later and dW a and dV a are the Wiener 
process defined to have following correlations: 

(dW a (x, t)dW b (x\ t)) = 25 ab 5 n (x - x')dt, 

(dV a (x,t)dV b (x',t)) = 25 ab 5 n (x - x')dt, 

(dW a ) = (dV a ) = (dW a dV b ) = 0. 

The probability distribution P(q,p,t) in the phase space is also defined through 
correlation functions as 



(•••) = J DqDp---P(q,p,t). 



Similar calculation to the case in the configuration space leads the 
Fokker-Planck equation 



dP 
~dt 



-dq A \dq B dq B J dpA^-dps dps 

where L is the Lagrangian in the phase space defined as 

L = ~ip A d q A + H. 

The Fokker-Planck equation has a solution in the equilibrium limit with which the 
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correlation function is rewritten as 



(•••) = J DqDp- • -exp(-L). 



The constant kernel C, introduced for dimensional consistency between the 
Langevin equations (2.6), has the same dimension as t -1 . It is fixed if we demand 
that the Langevin equations in the phase space (2.6) are reduced to the Langevin 
equation in the configuration space (2.2). For the purpose, we put the Hamiltonian 
(2.4) into Eqs.(2.6): 



dq A = -iG A *d Q p B 



AB t 



qAB qqCD 



BGcd 



+ G AB d l (G BC d l q C ) + 
dp A = C 



dQ AB. 



dt + E A dW a , 



dq B 

iG AB d q B - PA \ dt + VCE A a dV a . 



(2.8) 



We eliminate the momenta in Eqs.(2.8) and obtain the Langevin equation of q A 

E B a dV a dp B , G AB r dGcD o c. 



d q A = -iG AB d Q {iG BC d Q q u + 



c 



Vcdt cdt 



) 



-{- 



dq 



B 



-diq^diq 



D 



+ 



dG CD 



iGcEdQ(f + 



E , E c a dV a dp C 



dq B 

+ G AB d l {G BC d l q C ) + 



Vcdt cdt 

dG AB 



\cn a F , E D dV b dp D 
){tG DF d q + 



Vcdt cdt 



)} 



dq 1 



dt + E A dW a . 



(2.9) 



Taking the correlation of the Wiener process dV a in the limit dt — > 0, Eq.(2.9) 
should be identical with (2.2). It is realized if 

2 



C 



dt' 



(2.10) 



which is also supported in view of the dimension. We therefore conclude that the 
Langevin equations of the canonical variables q A and pa are 



dq A = -G AB ^dt + ^^dt + E A a dW a , 



dq 1 
dL 

dp A = -%Gabt\ h 

OPB 



dq B 
~2 



(2.11) 



E A a dV a . 



The drift term (the first term on the RHS) of the latter equation of (2.11) being 
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very powerful because of the order of dt, pa will reaches an equilibrium state in 
a very short time. It justifies the manipulation to connect Eq.(2.2) with Eq.(2.9); 
taking the correlation only of dV a accords to thinking the system at a finite time. 
The discussion on the constant kernel in this section may also be applied when we 
write the Langevin equation of an auxiliary field, the kinetic term of which does 
not exist in the Lagrangian. 



4. §3. FIRST CLASS CONSTRAINTS AND STOCHASTIC CONSIS- 
TENCY CONDITIONS 

In this section we consider the system, the Lagrangian of which is L(q I , doq 1 ) 
(I — , N), with FCC 



The conjugate momentum pj and the Hamiltonian H are defined through Eqs. 
(2.3) and (2.4) respectively and the metric of the system is assumed to be Sjj 
for simplicity. Following the canonical approach the Hamiltonian is generalized to 
include the FCC as 



where A a is a multiplier field, which cannot be determined without gauge fixing in 
the canonical formalism. We define the Langevin equations accordingly as 



F a (q,p)=0, 



(a = 



M). 



(3.1) 



H G = H + \ a F a , 



(3.2) 




where the singular Lagrangian Lq is 



Lg = -ip 1 dog 1 + H G . 



(3.3) 



It is convenient to define symplectic variables 



<f> l = { 



,q\---,q N }, 
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dY i = {dV 1 , • • • , dV N , dW 1 , • • • , dW N }, 
with which we can rewrite the Langevin equation in a simple form 

dtf^-^dt + dY*. (3.4) 

When we quantize gauge field with the Langevin equation in configuration space, 
we do not need to fix the gauge freedom; the Faddeev-Popov ghost effects are 
automatically carried into correlation functions and divergence, though appearing 

[3] 

in the correlation functions, does not in physical quantities. The Langevin equation 
(3.4), however, is not suitable for perturbation because of the multipliers, and 
neither is the Fokker-Planck formulation. The Fokker-Planck equation derived 
through the Langevin equation (3.4) is 



dP d fdL G , d 
~dt 



whose solution in the equilibrium limit is 

P = exp(-L G ). (3.6) 

We cannot obtain proper correlation functions either by using the distribution 
(3.6). 

We need to change the Langevin equation (3.4) into the form which fits for 
our handling. We, following the case in the configuration space , [4M try to impose on 
each FCC (3.1) the stochastic consistency condition 

dF a = (3.7) 

and hence determine the multiplier A a . Our strategy may seem strange since the 
multipliers for FCC cannot be fixed in the canonical formalism. We will however 

rgj 

accomplish the purpose as the following. If we take Ito's formula into considera- 
tion, the stochastic consistency condition (3.7) becomes 

= dF a 
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dF a • 1 d 2 F a ■ ■ 

dF" <(SL , dF b \ u j, ,1 

where {d4> l d^)o[dt) ls the correlation whose order is rft. It is not 28 %J dt since A a 
contains the Wiener process dY % and hence after setting the multipliers concretely 
we will compute it! 51[121 We change Eq. (3.8) into 

^ dF a dL 1 d 2 F a IJ±ij±js dF a dY i 



where 



<90 l .90* v ; 



Assuming that -D a& has the inverse D^, we obtain the definite form of A a 

1 r 

A a = D~ h \ — -—r —7 + 



_ l( dF b dL 1 <9 2 F b /JjLijJL ,- x aF 6 ciFM , oi1 . 

" ^ 1 " + 2^^w (rf0 d ^ )o <*) + ^^r/- (3 - n) 



The Langevin equaiton (3.4), into which Eq.(3.11) is put, becomes 
( f) T 1 1 f) P a cP 1 F b 

# = A " J {-a?* + dYi ) - 2W D: *WW^ W)o ^' (3 ' 12) 

where K iJ is the projection operator which extracts the quotient space composed 
of only the independent variables: 

K '■' = 5 ij /,"'•'. 

90* ab d(f>> 

In this stage we can compute {d(p' l d(j> l )o{dt) by using the Langevin equation (3.12) 
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and finally obtain the Langevin equation without the multipliers: 



* = ^{"a?* + dy 1 - W D ~ J WW K ' kdt ' (X13) 

which is the desirable form for calculation and of course satisfies Eq. (3.7). Im- 
posing the stochastic consistency condition thereby enables us to determine the 
multipliers for the FCC, which fact suggests that the stochastic consistency condi- 
tion may be equivalent to some gauge fixing condition. 

The Fokker-Planck equation, expected to bring us more information, is 

dP(M d ti fdL d*F* ,dF b d X 

dt ~ dft \d<j)j d(f>> d(f) k b d(f) k d(j)i J [ } 

Taking the FCC (3.1) into consideration, we write the equilibrium solution as 

P = J D\VDexp(-L- X a F a ), (3.15) 

D = det(D a6 ), 

where appears the determinant factor, which does not exist in the previous solu- 
tion (3.6). We hence expect that the FCC (3.1) and the gauge fixing conditions, 
which we will write as C a = 0, due to the stochastic consistency conditions (3.7), 
constitute the second class constraints whose Poisson brackets satisfy 

det{F a ,C b } pB = v^D. (3.16) 

The probability distribution (3.15), however, seems to include no gauge fixing term 
and diverge in the same manner as the previous one (3.6). Is the guess wrong ? 
If not, we must find the finite solution of the Fokker-Planck equation (3.14). The 
answer will be obtained in the next section. 
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5. §4. THE STOCHASTIC CONSISTENCY CONDITION AS GAUGE 
FIXING 

Our purpose in this section is to confirm the above mentioned conjecture by 
investigating in more detail the Langevin equation (3.13) and the Fokker-Planck 
equation (3.14), for which reason it is necessary to introduce new symplectic vari- 
ables and divide them into constraint and independent variables. We write the 
new variables as 



$m = { 



n 1 



n 



jV 



,Q\---,Q N }, (// = !,••-, 2N) 



(4.1) 



where Q and n are canonical conjugate of each other. Vielbein fields are, connect- 
ing the new variables with the old ones, defined as 



(4.2) 



E'E^ = 5£, E\E^ = 5, 



(4.3) 



Besides the definition, since Q and n are the canonical variables, the vielbein fields 
must satisfy the following relations: 



E 



a+N 



I+N 
E a+N 



E a ,i, 



J — —^a,I+N-> 

TP TP a 

&a+N,I+N — & I, 

TP - TP& 

&a+N,I - ~E I+N . 



(a 



N;I =!,•••, N) (4.4) 



[131 

Each of the constraints F a , assumed to be the first class, may be transfered to the 
new momentum I\. a / \/C, that is, we set the vielbein fields (4.2) with the left-hand 
indices /i — a (a — 1, • • • , M) as 



E a A 



dF a 



F -D^ dFb 



(4.5) 
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From Eqs. (4.4) and (4.5) we set those with fj, — a + N: 

_! dF b 

oh Q^I+N ' 

1 dF b 



E a+N j = —D 



rpa+N 
^ I+N 

E a +N,I 
E a +N,I+N 



ab -q^T' 

dF a 
dF a 



(4.6) 



There are some more relations which Eqs. (4.4) and (4.5) yield: 



& i-fr i — &a+N,i&b+N,i — u , 



E a+N E b+N = D -l 
i i ab ' 



E^gb+N _ E a jfib+N — E^Eb+Nj = E aj iEb + N,i = 0, 



R 

K, 



E iE a ,j, 
? a+N 



(4.7) 



ij ~ E a+I jE a+N j + E A E A j + E A+ i E A +N,j, 
(A = M + 1,---,N). 

The third equation of (4.7) is obtained with the help of the assumption that the 
constraints (3.1) are the first class. Moreover E A and E A+ f chosen orthogonal to 
each other and to the other vielbein fields, the metric of the manifold spanned by 
the independent variables $ a+iV and & A+N is 



( 9ab \ 
D ab 



V g AB ) 



{n'i/ = M + 1, • • • , 2N) 



(4.8) 



with 



g AB = E A E B { = EA+N,iE B +N,i, 

9AB = E Ai E Bji = E A+N lE B+N t . 
We thus obtain the Langevin equations of the new independent variables: 



+ 



K lJ dt 
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r F)T ■> f)F^ 

= G "'"'{-M* dt + E "'' iY '} + luf K " dt - (49) 



while the constraint variables of course satisfy d& a = 0. Here we assume 

V "' B "•' s 7BM? ( ^ E = °' (4 ' 10) 

where we use det(G^v) = D and V ^ is the covariant derivative in the manifold 
spanned by the independent variables, which assumption naturally leads metric 
condition V il >G pjV = 0. The Langevin equation (4.9) with Eqs. (4.8) and (4.10) 
is thereby reduced to 



which is GCT covariant in the manifold. 

The Fokker-Planck equation made up of the independent variables is 

d ,/r <9L 1 dVD d 



dP 
~dt 



which is no more than the rewriting of Eq. (3.14) and has only the same solu- 
tion. Calculating the expectation value only of physical quantities, we can however 
change the Fokker-Planck equation (4.12) into the desirable form which has a finite 
solution. Physical quantity must not depend on the value of the multipliers A a , 
that is, the Poisson bracket of FCC and the physical quantity / should be zero; 

{/> fl }p. B . = q§L = °> ( 4 - 13 ) 

which yields 

{df) = { d^ d + d^B 9 + dQA+NQB+NSAB). 

The Fokker-Planck equation (4.12) hence loses the gauge variables Q a+N and is 
reduced to 



dP 
~dt 



d AB [dL_ 1 dVD d 



} 
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d ( dL 1 dVD d -i 



(4.14) 



Finding the finite solution identical with corresponding path integral distribu- 
tion, we decompose the vielbein fields as 



% — e m %■> 



(/2 = M+l,...,2iV) (4.15) 



with 



J* f 



p A _ A+N _ a+N _ a+N _ « 



We then define new variables C a+N (a = 1, • • • , M) as 

dC a+N 



~QQa+N = e a+iV ' (4-16) 

with which we can write the equilibrium solution of the Fokker-Planck equation 
(4.14), setting C~ a+N = 0, as 

P = J D\D\VDexp(-L + \ a F a + \a +N C a+N ). (4.17) 
Taking notice that 

det{C &+N , F a }p. B . = deb(eJ» N ) = Vfi, 
the Fokker-Planck distribution (4.17) is identical with the path integral distribution 
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which we obtain by imposing the gauge fixing conditions 

C a+N = 0, (a = l,---,M) 

upon the FCC 

F a = 0. 

We thereby insist that our conjecture in §3 has been confirmed. 

6. §5. CONCLUSION 

In this paper we have formulated the SQM-ps on the classical basis of the 
canonical formalism; the Langevin equations in the phase space reduce to Hamil- 
ton equations in the classical limit. Furthermore our prescription of introducing 
constraints into the Langevin equations is consistent with Dirac's formalism. The 
Langevin equations with FCC, however, do not need to contain ordinary gauge 
fixing conditions, in place of which stochastic consistency conditions have been 
imposed and enabled to determine the Lagrange multipliers for the FCC. It has 
suggested that the stochastic consistency conditions should correspond to some 
gauge fixing. More information has been brought through the equilibrium Fokker- 
Planck distribution: After imposing the stochastic consistency conditions upon 
the FCC, we can obtain the finite values of physical quantities though the other 
quantities may diverge in the limit t — > oo. We therefore conclude that impos- 
ing stochastic consistency conditions on FCC, being a little differnt from ordinary 
gauge fixing procedure, makes the FCC change into second class constraints. 
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8. APPENDIX 



Here we show that our formulation is applicable to the systems with second 
class constraints. Following our prescription with gauge-fixed Hamiltonian 

H G = H- \ a F a - ~\~ a+N C~ a+N 
= H- A a F\ 
(a=l,..-,2M) 



we can straightforwardly reach the Langevin equation without the multipliers: 



where 



d<p l — K \ —^Yjdt + ^ > - —D &h _K* dt, 



-ab «9F a <9F b 



,90* ab <90J " 

We can easily set up the Fokker-Planck equation, too: 

dt <90* \d(j)3 d^d<t> k ab d<t> k oft) w ; ' 

which has an equilibrium solution 

P = J D\VT)exp(-L + A a F a ). (Al) 

The solution (A.l) should be identical with corresponding path integral distribu- 
tion, which we can confirm with the following equations: 

d*{C*»**} M .=d*(^) 
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